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Abstract 

In this paper analytic formulas for electricity derivatives are calculated. To this end, we assume that 
electricity spot prices follow a 3-regime Markov regime-switching model with independent spikes and drops 
and periodic transition matrix. Since the classical derivatives pricing methodology cannot be used in case of 
non-storablc commodities, we employ the concept of the risk premium. The obtained theoretical results are 
then used for the European Energy Exchange (EEX) market data. The 3-rcgime model is calibrated to the 
spot electricity prices. Next, the risk premium is derived and used to calculate prices of European options 
written on spot, as well as, forward prices. 
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1. Introduction 

Deregulation of electricity markets has led to a substantial increase in risk borne by market participants. 
The often unexpected, extreme spot price changes range even two orders of magnitude and can cause severe 
financial problems to the utilities that buy electricity in the wholesale market and deliver it to consumers 
at fixed prices. The utilities and other power market companies need to hedge against this price risk. A 
straightforward way to do it is to use derivatives, like forwards and options. It is exactly the aim of this 
paper to price commonly traded electricity derivatives. 

Before calculating a price of a derivative, a proper model for the underlying asset has to be chosen. There 
are two approaches comm on for the electricity ma r ket. The first one is to start wi t h spec i fying the forward 
price dynamics (see e.g. Idewow and Strickland Il999t iBenth and Koekebakkerl l2008ct iBierksund et all 
2010). Such approach is useful if only derivatives written on forwards are to be considered and a link 



between forward and spot prices is not important for modeling issues. Obviously, the spot price can always 
be derived using the fact that the forward and spot prices should coincide at the forward settlement. 
However, the complexity of the spot price dynamic s is then usually neglected. The second approach is based 
on defining the spot pr i ce dy namics first (see e.g. iLucia and Schwartz! . l2002[ iMiltersenl l2003t IBenth et all 
2003t iBierbrauer et all 120071 ) and then modeling a link between the forward and spot markets. U sually, a 
convenience yield or risk premium notion is used ( Benth et all l2008at iGemanl . l2005t IWeronl . 120061 ) . Using 
such an approach allows to price derivatives written both on the spot and the forward price. Moreover, 
a relation between spot and forward prices is taken into account and the lack of forward price data is no 
more a limitation. Here, we use the latter approach and describe the spot price dynamic by a Markov 
regime-switc hing (MRS) model with ind ependent spikes and drops and periodic transition matrix that was 
proposed by Ijanczura and Weronl (J2010I) . 

After specifying a model we have to choose derivatives pricing methodology Classical approach used in 
financial and commodities markets is based on the no-arbitrage assumpti on and construction of a strateg y 
replicating a future payoff (or equivalcntly finding a martingale measure, iMusiela and Rutkowskil (J1997I )). 
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However, such approach fails in case of electricity due to very limited storage possibilities. Therefore, instead 
of using a martingale approach, we employ a concept of the risk premium/market price of risk and find such 
pricing measure that yields the observed forward market prices. With such methodology we are able to 
derive forward prices from the spot price model and also to find explicit formulas for premiums of European 
options written on spot, as well as, on forward prices. 

The paper is structured as follows. In Section [5] we introduce the Markov regime-switching model used 
for electricity spot price dynamics. Next, in Section [3] we explain why the classical derivatives pricing 
approach fails in case of electricity derivatives and, as a solution, we describe the 'risk premium' approach 
in case of the considered model. The obtained results are then used to derive analytical formulas for prices 
of electricity derivatives in Section [4] Finally, in Section \E\ we use the obtained theoretical results to price 
derivatives using the European Energy Exchange (EEX) data and in Section [5] we conclude. 



2. The model 

Let us first recall the main stylized facts about electricity prices. Electricity prices highly depend on the 
actual demand/consumption. Obviously, the latter is varying during a year due to the changing weather 
conditions and throughout the week or day due to the business cycle. The same long-term (yearly) and 
short-term (weekly/daily) seasonality is recorded for electricity prices. The second apparent feature is the 
very high volatility of electricity prices and unexpected, usually transient, dramatic price changes called 
spikes or jumps. The spot price may rise for a few hours to even two orders of magnitude of the standard 
prices and then fall back to the normal level. What makes electricity market completely different from 
other financial or commodities markets is that the electricity prices may, as well, abruptly fall down yielding 
negative values. Finally, electricity prices are mean-reverting, meaning that in long time period they move 
back to some equilibrium level. 

Seasonali ty is usually removed from the analyzed prices prior to mo deling by fitting some p eriodic func- 

tion like sine (Pilipovid.ll998tlCartea and Figueroal 120051 ; iDe Jona 120061) or piecewise constants ( Lucia and Schwartz , 
2002[ iKnittel and Roberts!. 120051 ). Alternatively, some smoothing technique like wavelets or moving average 
can be used ( Weronl . 120091 ) . The mean-re verting proper ty is typically modeled with some mean- reverting 
processes, like e.g. AR(1) time series or the lVasicekl ( 1977 ) model. The most challenging for modeling, and at 
the same time the most important for risk management, are the pri ce spikes. One approach is to incorporate 



a jump comp onent into a standard mean- reverting diffusion model (jDengl . Il998t ICartea and Figueroal 12005 



Weronl . 120081 ). However, in the resulting jump-diffusion models there is a problem of how the price after a 
spike revert back to the normal level. Nor an immediate negative jump, nor mean-reversion pulling back 
the prices to the normal level, yields a flexible tool for modeling consecutive spikes. Another possibility are 
the Markov regime-switching (MRS) models in which the prices might stay in the excited (spike) regime 
with some probability. Hence, MRS models allow for modeling consecutive spikes in a very natural way 
and seem to be a reasonable choice for electricity price dynamics . To our best knowledge MRS models 
were first applied to electricity prices bv lEthier and Mountl ( 1998T ) who used an AR time series with pa- 
ra meters depending on the actu al re gime. A MRS m odel with independent spikes was later introduced 
by iHuisman and De Jonj (120031) and De Jong J 2006 \. N umerous attempts i mproving statistical proper 



ties of the m odel (see e.g. IHuisman and Mahieul . 12003c I Mount et al. 



Weron . 20091) or includ i ng so me exogenous factors (see e.g. JHuismaiu 



r. 



20061: lHaldrup and Nielse: 
20081 : Kanamura and Ohashil 



,:: 



2006 



2008 



Karakatsani and Bunnl . 120081) were late r proposed. Here, we focus on a 3-regime model with independent 
spikes and drops introduced recently bv ljanczura and Weronl ( 2010I ). 

Having in mind the above mentioned features of electricity, we let the electricity spot price be given by 



P, 



9t 



X, 



(1) 



where gt is a deterministic seasonal component and X t follows the 3-regime model with independent spikes 
and drops. Namely, 



X t = 
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X 
X 
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where b denotes the base regime (describing the 'normal' prices) , s the spike regime (representing the sudden 
upward price jumps), while d stands for the drop regime (responsible for the sudden price drops). Further. 
\t\ denotes the integer part of t, Rk is a discrete-time Markov chain defined by a time-varying (periodic) 
transition matrix P(fc) 

(Pbb(k) p bs {k) p bd {k) \ 
p ab {k) Pss {k) Psd (k) , (3) 

Pdb(k) p ds {k) Pdd(k) J 

for k € {0, 1, 2, ...}. The base regime dynamics is given by the IVasicekl (119771 ) model: 

dX t , h = (a- pX t , b )dt + a b dW t , (4) 

having unique mean-reverting solution of the form: 

X t = Xoe-* + % (1 - e-*) + a f e -^^dW u , (5) 

P Jo 

where Wt is a Wiener process (or Brownian motion), (3 is the speed of mean-reversion, % is the long-time 
equilibrium level and a is the volatility. The spike regime values (Xo tS , Xi iS , -X^ lS , ...) constitute an i.i.d. 
sample from the c s -shifted log-normal distribution, i.e. : 

]n(X M -c,)~JVG* aj c*), X kiS >c s (6) 

while the drop regime values (^o,d,^i,d,^2.d, •••) form an i.i.d. sample from the inverted c^-shifted log- 
normal distribution defined as: 

ki(-X M + c d ) ~ N(fi d , al), X k4 < c d - (7) 

Observe, that in the model defined by ^ the price process can jump to a different regime only at discrete 
time points t = 0, 1, 2, .... This is motivated by the fact, that, even though the price is a result of continuous 
bidding, electricity spot price is typically settled for contracts with some delivery period, usually an hour. 
Hence, the change in electricity spot price dynamics may occur only in discrete time points. 

Finally, let (Q, J 7 , Q) be a probability space with filtration T t generated by the processes X t and Ry t \ 
and assume a constant continuously compounded interest rate r. 

3. The risk premium 

The classical option pricing approach is based on the no-arbitrage assumption, which implies that the fair 
price of a derivative is the discounted expected future payoff under a martingale measure ( Harrison and Pliskal 



Il983l ). If the market is complete, i.e. any contingent claim can be replicated with a self- financing strategy (it 
is attainable) , there exists a unique martingale measure. However, due to the non-storability of electricity 
a derivative written on the spot electricity price cannot be replicated with a portfolio consisting of the un- 
derlying instrument and a financing debt account. As a consequence, the market is incomplete. Moreover, 
since electricity cannot be traded in the usu al way (once purcha sed has to be consumed), the only tradable 



asset in the spot market is the bank account (IBenth et al\ , 120031) . Recall that here we assume a continuously 



compounded constant interest rate r. Obviously, the discounted value of the bank account is a martingale 
under any measure equivalent to the actual (also called the objective or statistical) measure Q. Therefore, 
the market is arbitrage-free but there is no unique martingale measure. An additional criterion has to be 
used in order to select a pricing measure. 



H ere we use an approach based on the concept of the risk premium (see e.g. IBenth et all , 12008a: G eman 



2005t IWeronl 120061) , which is defined as a reward for investing into a risky asset instead of a risk- free one. 



Note, that a related notion is that of the market price of risk, which can be seen as a drift adjustment in 
the dynamics of an a sset to reflect how investors are compensated for bearing risk when holding the asset 
( Benth et all [2008bJ). The idea of the 'risk premium' approach is to choose a martingale measure that is 



consisten t with the prices of forw ard contracts quoted in the market. A similar approach is used in the 
weather (jBenth and Benthl . 120071 ) or interest rate derivatives context and is based on calibrating the model 
to the initial yield curve (see e.g. iHull and White! . Il993t lBi0rkl . ll997l ). 

Recall, that the arbitrage-free price of a forward contract should be equal to the expected future spot 
price under a pricing measure, namely 

fl = E\P t \F ), (8) 

where /q is the price at time of a forward contract with a delivery at time t, Pt is the electricity spot price 
and E x (-) is the expected value with respect to the pricing probability measure Q x , equivalent to the actual 
measure Q. In an incomplete market relation ([8]) does not yield a unique forward price, as it is dependent 
on the choice of Q x . Here we choose Q x such that relation ([5]) is consistent with market data, i.e. E x {P t \J : tf) 
is calibrated to the quotations of forward contracts. In other words, we choose the pricing measure that is 
used by the market. In the following we will assume that the measure Q x is the probability measure under 
which the drift of the base regime process is parametrized by a function X(t) (i.e. the market price of risk) 
chosen so that E x (P t \J r o) yields the market forward price /q. 

Before we find the measure Q x , we give a brief explanation of how to calculate the risk premium in the 
3-rcgimc model dU-©. Assume, that the forward price /q is given for any maturity t and the spot price 
model parameters 9 = (a, j3, ab, /U S) o~ s , c s , \id, &d, Cd, P) are known. 

The risk premium RP is defined as 



RP(t) = E(P t \ Jb) - fl 
where /q is the market price at time of a forward contract with delivery at time t. 



(9) 



Remark 3.1 It should be noted that some authors (see e.g. \Evdeland and Woluniea . \2003 ) define the risk 
premium as the difference between the forward price and the expected spot price, i.e. —RP(t). 



,(*) 



Let p\- = P(Ry t \ = j\Ro = i) denote the probability of switching from state i at time to state j at time 
\t\ . For a constant probability matrix P it is given by the ijth. element of the \t\ th power of the transition 



matrix, i.e. p\- = (pL*Jj ... For a time- varying probability matrix it is given by p\. = ( Ilfc=o -^CO 

In order to simplify the derivation, in the following we assume that P(Rq = b) = 1 or equivalently 
Xq = Ao.b a.s., i.e. at time the process X t is in the base regime with probability 1. 
Now, we can derive a formula for the risk premium. Observe that 

E{X t \T ) = P(R [ti = b\R = b)E(X tlb \? ) + P(R vn = s\R Q = b)E(X ltls \T ) + 

+P(R [ti =d\R = b)E(X [tl jT ). 

Recall, that -^|tj, s and -X"|t|,d arc random variables independent of J-q. Hence, E(X^ t jj\J : 'o) = E(X^j) for 
j = s, d. Moreover, from the assumption of Xq = Xo.b we have that E{X t ,b\Fo) = E(X tt b\Xo.b). Hence 



(*); 



E(X t \F )= Pb l ! E(X tt b\X 0t b) +pH'E(X ltU ) +p^E(X ltitd ). 



(10) 



As a consequence, from ([T]), (|SJ| and (|10p . the risk premium in the 3-regimc MRS model defined by equations 
©-©is given by: 



RP(t)=P% 



x e 



-pt 



a 



(1-e-*) 



■^{^ 



+ Cs)+ Pbd -e w+5 ' d +Cd)+9t- /o 



where xq is the stochastic part of the price observed at time and /q is the market forward price. 



Remark 3.2 Observe that, if assumption that Xq = Xq^ is not satisfied we have: 

oo 

E(X tyb \T Q ) = I{ Ro=b }E(X t ,b\Xo.b) + 2^ ^{R ^b,R-i^b,R^b,...,R- k+1 ^b,R- k =b}E(X t , b \X_ k+1 . b ), (11) 



fe=l 



where a negative time index is used for the historical (i.e. before the moment of valuation t = 0) values of 
the process and for u < t E(X t , b \X u , b ) = X^e - ^* - ") + § (l - e'^"^). Note, that formula ^ is a 
consequence of the fact that the base regime values become latent if a spike or drop occurs. Moreover, 



P(R lti =j\T )= Y, hR = l} P(R[ti = j\Ro 

ie{b,s,d} 



(/) 



i )= J2 l {Ro=i}Pij 

ie{b,s,d} 



Thus, the risk premium calculation and all of the following results can be generalized to the case Rq =/= b. 



4. Electricity derivatives pricing 

4-1. Options written on the electricity spot price 

Now, we turn to pricing of a European call option written on the electricity spot price. Recall, that 
the European option is a contract that gives the buyer the right to buy /sell the underlying commodity at 
some future date t (called maturi ty) at a certain price K (called the strike price). First, we find the pricing 
measure Q x . Like iMertonl (119761 ) in the context of jump-diffusion processes we assume that the dynamics 
of spikes and drops are the same in the actual and pricing measures. We start with finding the spot price 
dynamics under A parametrization. 

Let \(u) be a deterministic function square-integrable on u € [0,T max ], where T max is a time horizon 
long enough to contain all maturities of derivatives quoted in the market, and introduce a new process W x : 



W t x = W t 



AH 

a b 



du, 



(12) 



where a b is the volatility of the base regime. From the Girsanov theorem we have that W x is a Wiener 
process under a new measure Q x defined as 

dQ x 



dQ 



exp 



Ob 



dW„ 



AH 
<Jb 



du 



with the filtration J^ , being the natural filtration of the process Wt- 
Now, the base regime process X t .b can be rewritten as: 



dX tib = [a- X(t) - pX tlb ]dt + ObdWl 



(13) 



(14) 



and the expected future spot price is given by: 



(') 



E\P t \F ) = p% 



Xne 



-0t 



+*2 (- 



C<I 



-0t\ 



9t- 



(t-u) 



X(u)du 



.^(e^M 



C S + 



(15) 



The function A(i) can be calibrated to the market forward prices so that E x (P t \J : o) = /q, e.g. by using 
some fitting procedure (like the least squares minimization). Alternatively, one can find the risk premium 
and then use the relation between the market price of risk X(t) and the risk premium: 



/o 



fi{ - t - u) \{u)du = RP{t), 



(16) 



which is a simple consequence of the fact that RP(t) = E(P t |.Fo) — E x (P t \Po) , formula ([±3]) and Ito's lemma. 
Now, the price of a European call option written on the electricity spot price can be derived. 



Option price formula. If the electricity spot price P t is given by the MRS model defined by equations (HJ-G]), 
then the price of a European call option written on P t with strike price K and maturity T is equal to: 



where 



C T (K) 



-rT 



P^CtAV+P^CtAV+p^CtAK) 



C T . b {K) 



cxp 



(K' - mf 
2s~ 2 



+ (m-K') 



1-$ 



K' - m 



(17) 



(18) 



CtA k ) = I{A"> Cs }< exp/x s + ^ 



1-$ 



/ log(K'-c a )-ti s -cr 2 a 



-{K 1 - c s ) [1 - Fln^.^K' - c,)] \ + I { k><c 



exp /i 



cs - A" 



and 



CtA K ) = hK'<c d }' 



fM + Y' 



log(c d - K') -ii d -o-j 



<7,l 



+ (c d - K')F LN(lid(T 2 d) {c d - A') 



Further, K' = K 



Jo e 



/3(T-u) 



X(u)du, 



1 



-2^T\ 



and 
2 



9T , m = X e-? T + f (1 - e-^ T 7 
Alw(u,o- 2 ) is the cumulative distribution function of the log-normal distribution with parameters /z and a 

Note that, in order to make the exposition of the paper clear, the price derivation is moved to the 
Appendix. 

Here, we assume that the option is settled in an infinitesimal period of time [T,T + A]. However, in 
practice, the electricity spot price usually corresponds to a delivery during some period of time (e.g. an 
hour, a day) and, hence, the maturity of the option should be specified on the same time-scale. On the 
other hand, the analyzed spot price quotations usually represent some delivery period. For instance, if the 
considered data is quoted daily, as it will be in the empirical example of Section [5j then the maturity of the 
option would be also given in daily time-scale and would correspond to daily delivery. 



4-2. Electricity forward contracts 

Probably, the most popular electricity derivatives are the forward contracts. Recall that a forward 
contract is an agreement to buy (sell) a certain amount of the underlying (here MWh of electricity) at a 
specified future date. Settlement of the contract can be specified in two ways: with physical delivery of 
electricity or with only financial clearing. Both types of settlement are in the following called delivery. 
Denote the price at time t of a forward contract with a delivery at time T by / 4 T . Since the cost of entering 
a forward contract is equal to zero, the expected future payoff under the pricing measure should fulfill: 



what implies that 



E x (P T -fJ\F t ) = {), 



ff = E\P T \T t ). 



(19) 



(20) 



Observe, that now we define the price of a forward contract at any future date t. This is motivated by the 
fact that the valuation at time of an option written on a forward contract requires the knowledge about 
the forward price dynamics at the option's maturity t. 



Forward price formula. If the electricity spot price Pt is given by the MRS model defined by equations 
CO)-®, then the price at time t of a forward contract written on P t with a delivery at time T is given by 
the following formula 



i\ 



P(R m = b\T t ) 



E x {X t , b \T t )e- 



(T-t) 



(T-t\ 



P 



-P(R\t] = s|Ji)(e" 



P{R V t\ = d\F t ){c d 



,W+5»j 



e-"( T -")A(«)d« 



) + 9T, 



(21) 



where P{Ry T ^ = i\Ft) = Yl 



je{b,s,d} 



P(R m =i\R lti =j)l {R[t . =j} 



Note that in the above formula E (X t .b\J~t) is used, since this expectation depends on the state process 



value at time t. Namely, if Rt 



bthenE x (X t , b \T t ) 
E x (X t . b \T t 



X t .b = X t . On the other hand, if at time t a spike or 
_i) and again this expectation is dependent on Rt-i value. 



a drop occurred then E {Xtj,\Ft) 

A general formula for E x (X tjb \J r t) can be found using the same derivations as in Remark 13.21 

When deriving the forward price dynamics, we have to remember that the properties of the obtained 
model should comply with the observed market prices. One of the most pronounced features of the mar- 
ket forward prices is the observed term structure of volatility, called the Samuelson effect. Precisely, 
the volatility of the forward prices is quite law for distant delivery periods, however, it increases rapidly 
with approaching maturity of the contracts. Here, the forward price volatility is described by the part 
P{R[T\ = b\F t )E x {X tfi \J : t)e~ l5i - T ~ t *> of formula l[2"T]), Hence, it is specified by the volatility of the spot 
price base regime scaled with e _ ^ T_t ' and the corresponding probability of switching to the base regime. 
Observe that the scaling factor e - ^ T-t ) exhibits the Samuelson effect as it increases to 1 with t approaching 
maturity time T. Moreover the forward price volatility, again due to the scaling factor, is lower than the 
spot price volatility. This is in compliance with the behavior of the market spot and forward prices. 

Electricity forward contracts listed on energy exchanges are usually settled during a certain period of 
time (a week, a month, a year, etc.). Denote the price at time t of a forward contract settled during the 
period [Ti,T 2 ] by f\ u 2 '. Obviously, the latter is the mean price of forward contracts with delivery during 
the period [Ti,T2], namely: 



/, 



[Ti,T 2 



21 = / ' w(T 1 ,T 2 ,T)f?dT = [ 2 w(T u T 2 ,T)E x (P T \T t )dT. 



(22) 



where w{Ti, T 2 , T) is the weight function representing the time value of money. The form of w depends on 



the contract specification. For contracts settled at maturity we have w{T\, T 2 , T) 



T 2 ~T! 



while for instant 



settlement w(Ti,T 2 ,T) — -Ji\_ - t t 2 , where r > is the interest rate (JBenth et all . l2008ah . The price 
fl 1 ' can be obtained from formulas (|2~Tj) and (|22|) . Indeed, we have: 



f [TuT 2 ] = E x {XtblTt) f 2 w ( Tu T 2 ,T)P(R lT] =b\T t )e-^ T -^dT + 



f 2 U ;(r 1 ,r 2 ,r)p(i? m =b\T t ) 



- ( e-^> 



e-^ T - u ^X(u)du 



dT 



+(e^ 



'•+c t ) / w(T l7 T 2 ,T)P(R m = s\F t )dT + 
+(c d -e^+^) f 2 w(T 1 ,T 2 ,T)P(R lTl = d\F t )dT + f ' w(T u T 2 ,T)g T dT. (23) 



4-3. Options written on electricity forward contracts 

Finally, we find an explicit formula for a European call option written on a forward contract delivering 
electricity during a specified period of time. Observe, that the forward price f\ u depends on the spot 
price at time t and, as a consequence, also on the state process value at time t. 



We consider an option written on an electricity forward contract with settlement during a specified period 
of time, as it is the most popular specification of electricity options on energy exchanges. For example, in the 
EEX market there are options written on forward contracts with monthly, quarterly and yearly settlement 
periods. The maturity of such options is set to the fourth business day before the beginning of the underlying 
contract's settlement period. The EEX example will be examined in Section [5] 

Price formula for an option written on a forward contract. If the electricity spot price Pt is given by the 
model defined by equations (U}-([7]), then the price of a European call option with strike price K and maturity 
t written on a forward contract with delivery during the period [Ti,T 2 ] is equal to: 



where 
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w 



(T 1 ,T 2 ,T)P(i? LTJ =b\R\t } = fyer^-^dT, 



(24) 

(25) 
(26) 
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,-T, 



-P{T-u) 
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LtJ-fe+i 
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_ e w+M) i w (T 1 ,T a ,T)P(R m =d\R lti =i)dT+ I w{T x ,T 2 ,T)g T dT, (27) 
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+(e 



-p(T-u) 



X(u)du 



dT 



+ (c d 



1 2 
o^d+2 a d 



T 2 

) I w(T u T 2 ,T)P(R lTi = d\R lti = b)dT 

Ti 



Ti 



w{T l7 T 2 ,T)g T dT. (28) 



T x 



and Ct t b(K) is the 'base regime part' of the price of a European call option written on the electricity spot 
price with maturity t and strike K, see equation (|18[) with T = t. 



5. EEX market example 

Theoretical results from the previous Sections allow us to price energy derivatives. We assume that the 
electricity spot price follows the model specified by equations (HJ)-© with a periodic transition matrix and 



Table 1: Coefficients of the function Lt, see equation I I29I I, fitted to the EEX spot prices. 
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Figure 1: Left panel: EEX spot prices and the fitted long term seasonal component (blue solid line). Right panel: weekly 
periodicity calculated for the EEX spot prices. Note, that 1 denotes Monday, 2 - Tuesday, ..., 7 - Sunday, while the eight day 
of the week is used for the German national holidays. 



c s , Cd being the first and third quartile of the dataset, respectively. We use mean daily EEX spot prices 
from the period January 2, 2006 - January 2, 2011 (5 years and 261 whole weeks). In order to calibrate the 
model, we first remove the seasonal component. 

We assume that the deterministic function g t is composed of two parts: a long term trend L t and a 
weekly seasonality St- Since the valuation of derivatives requires forecasting the seasonal component, we 
model the long term trend by a sum of sine functions: 



Lt = (ai + a,2t) sin [2-7T (t + a^)} + (04 + a^t) sin [27rag (t + 07)] + as + cigt + a\ot 



(29) 



where t is in yearly time scale. Note, that the first component of the above sum is responsible for the yearly 
periodicity, while the second one captures seasonalities of different period than one year (here, we obtain 
nearly half- year period, see o,q in Table [1]). The function L t is fitted to the EEX spot prices using the least 
squares method. The obtained curve is plotted in Figure [TJ while the obtained a, coefficients are given in 
Tabled] 

The estimated long term trend is subtracted from the analyzed time series. Next, the short term 
seasonal component is estimated using the 'average week' method. Namely, we calculate the mean of prices 
corresponding to each day of the week (German national holidays are treated as the eight day of the week). 
The obtained weekly pattern is plotted in Figure [1] 

Finally, the deseasonalized prices are obtained by subtracting the long and short term trend from the 
EEX spot prices. Moreover, the data is shifted so that the minimum of the deseasonalized and the original 
prices is the same. The resulting time series can be seen in Figure [5] 

After removing seasonality, we are left with modeling the stochastic part X t . Here, we calibrate the 
3-rcgimc MRS model, see equations ©-([7]), to the deseasonalized EEX prices. To this end, we use a version 
of the Expec t ation- Maximization algorithm o f iDempster et al.l (|l977j ). which was applied to MRS models 
by iHamiltonl ( 19901 ) and was later refined by iKiml (J1994I ) . For the detailed description of the algorithm 
in case of the 3- regime model considered in this paper see the recent work of iJanczura and Weronl (J2012I ). 
The obtained parameters are given in Table O Observe high probabilities of staying in the same regime, 
ranging from 0.40 for the drop regime up to 0.97 for the base regime. Hence, the assumed model allows for 
modeling consecutive spikes or drops in a very natural way. The calibration results are plotted in Figure [2 
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Figure 2: Calibration results of the 3-regimc MRS model fitted to the EEX descasonalized prices. The prices classified to the 
spike regime i.e. with P(S) = P(Rt = s\Xi,X2, ...,Xx) > 0.5 are denoted by red dots, while the pricess classified to the drop 
regime i.e. with P(D) = P(Rt = d\Xi, X2, ...,Xt) > 0.5 are denoted with black x's. Additionally, in the lower panels the 
corresponding probabilities arc plotted. The estimated unconditional probability of spike occurrence is displayed in the bottom 
panel. Observe, that the highest probability of spike is obtained for the Autumn/Winter period, while the lowest for Spring. 



where additionally the regimes classification is illustrated. As we may observe spikes/drops occur usually as 
a series of high/low prices rather than separate outstanding observations. What is interesting to note, is the 
clear seasonal pattern in the estimated probability of spike, see the bottom panel in Figure O Indeed, the 
highest spike occurrence probability is obtained for the Autumn/Winter season, while the lowest for Spring. 
In order to validate the used MRS model, we apply a Kolmogorov-Smirnov goodness-of-fit test for 
the marginal distribution of the individual regimes, as well as, for the whole model. We use two testing 
procedures. The first one (called ewedf - equally weighted empirical distribution function) is based on 
classifying observations to the most probable regime, i.e. assuming that R t = i ii P(Rt = i\Xi, X2, ..., Xt) > 
0.5. As a consequence, the standard Kolmogorov-Smirnov goodness-of-fit test can be applied. The second 
one (called wedf ) utilizes a notion of the weighted empirical distribution function, where i-th observation is 
taken into account with weight proportional to the probab ility P(Rt = i\Xi,X2, ...,Xt)- For the detailed 
testing procedure derivation see ljanczura and Weronl(J2012l) . The obtained test p- values are given in Tabled 



Table 2: Parameter estimates of the 3-regime MRS model fitted to the deseasonalized EEX prices. 
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Table 3: p-values of the goodness-of-fit tests for the 3-regime MRS model fitted to the deseasonalized EEX prices. The test 
results for the ewedf, as well as, the wedf approach are provided. 
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Table 4: Specification of the monthly forward contracts listed on the EEX market on January 3, 2011. 



Name 


Settlement Price 


Ti 


T 2 


Feb- 11 


54.35 


1.2.2011 


28.2.2011 


Mar- 11 


51.64 


1.3.2011 


31.3.2011 


Apr- 11 


48.07 


1.4.2011 


30.4.2011 


May- 11 


45.53 


1.5.2011 


31.5.2011 


Jun-11 


48.50 


1.6.2011 


30.6.2011 


Jul- 11 


49.00 


1.7.2011 


31.7.2011 



Recall, that p- value higher than 5% means that we cannot reject, at the 5% significance level, the hypothesis 
that the analyzed dataset was driven by the assumed model. As all of the obtained p- values are higher than 
5%, we cannot reject the considered 3-regime MRS model as a proper one for the analyzed dataset. 

Before we start with the valuation of derivatives we have to find the risk premium and the function A, 
see equation (TIB)) . We use monthly forward contracts listed on the EEX market on January 3, 2011, i.e. on 
the day directly following the calibration period. The prices, as well as, the delivery periods of the analyzed 
forward contracts are given in Table |4] 

Assume that X(t) = Xit + A2. Since we analyze monthly contracts, the risk premium should be also 
calculated on the monthly basis, i.e. instead of equation Q we use: 

RP{T U T 2 ) = I £ E ( P ^) - /F 1,T2] . (30) 

Note that the average monthly expected spot price is calculated as an arithmetic mean instead of an integral, 
because the analyzed spot prices are quoted is discrete time (on a daily basis). For the same reason Ai and 
A2 are found by fitting: 

Ai ££ Tl pff(* - (1 - e~*)//3) + A 2 Ef: Tl pff(l ~ e-*) 

m _ Ti + i) - RP(T U T 2 ). (31) 

See equation (|16[) for the comparison with the continuous time scale. The values of the risk premium obtained 
from contracts with different delivery periods are plotted in Figure[3](for contract specifications see TableU]). 
Observe a strong evidence for the negative risk premium, especially for contracts with approaching delivery 
period. For the contracts with more distant delivery period the obtained risk premium is less significant. 
Using the least-squares minimization scheme we get A = 0.0084i— 1.8387, see the red dashed line in Figure 
[3] for the plot of the function fitted to the risk premium (i.e. the left hand side of equation (|3l"j) ). 

Now, we can derive the price of a European call option written on the electricity spot price. Assume 
that the interest rate r is equal to 0. The option prices obtained in Section |4. II with different maturities t 
and strike prices K are plotted in Figure [5] Obviously, the lower is the strike price, the higher is the call 
option price. What is interesting to note, is how the option price depends on the maturity tenor t. Observe 
the clear seasonal pattern of option prices, both on the weekly and the long-term level. The long-term 
seasonality is caused not only by the deterministic component g t but also by the periodic transition matrix 
allowing for varying spike (drop) probabilities during the year. Recall that in the EEX market the spike 
probability is the highest in Autumn/Winter and the lowest in Spring, see Figure [21 

Now, we derive the prices of European call options written on monthly forward contracts. Using the 
results obtained in Section 14.31 we calculate the price of an option written on the forward contract with 
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Figure 3: The risk premium obtained from monthly forward contracts with different delivery periods (blue lines), as well as, 
the fitted function (see equation II31I plotted with red dashed line. 
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Figure 4: Prices of a European call option written on the electricity spot price for different maturities T and strike prices K. 
The prices are calculated on January 3, 2011. 
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Figure 5: Prices of a European call option written on the electricity forward contract with settlement in February 2011 for 
different maturities t and strike prices K. The prices are calculated on January 3, 2011. 

Table 5: Prices of European call options written on monthly forward contracts with strike price K = 30. The prices are 
calculated on January 3, 2011. 



Month of delivery [2\, T2] Options maturity t 



Options price Cfl Tl ' T2l (K) 



February 

March 

April 

May 

June 

July 



26.1.2011 
23.2.2011 
28.3.2011 
26.4.2011 
28.5.2011 
27.6.2011 



13.5449 
9.1074 
4.2258 
2.6826 
3.6148 
7.1913 



settlement in February 2011. The results are plotted in Figure [2] In order to check how the option price 
varies according to the delivery period, we calculate the prices of options written on monthly forward 
contracts with deliveries within the next 6 months (i.e. February 2011 till July 2011). According to the 
products specification in the EEX market, the maturity of the options is set to the fourth business day 
before the beginning of the delivery period. The obtained option prices are given in Table [5] Similarly, as 
in the case of options written on the spot price, we observe the lowest option prices for settlement during 
the Spring months. 



6. Conclusions 

In this paper we have derived premiums of European options written on electricity spot, as well as, 
forward prices. We assumed that electricity spot prices can be described by a 3-regime MRS model with in - 
dependent spikes and drops and periodic transition matrix, proposed earlier bv lJanczura and Weronl ( 2010I ). 
The forward prices were then derived using the risk premium approach and fitting the model-based prices 
to the observed forward curve. Next, using the spot and forward price dynamics we calculated prices of the 
corresponding European options. 
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The assumed model was then calibrated to the spot prices from the European Energy Exchange. We 
have validated the model choice by performing a statistical goodncss-of-fit test. Next, using monthly forward 
contracts quotations we have calculated the risk premium. We have obtained negative values, especially 
significant for contracts with approaching maturity. For contracts with distant settlement the risk premium 
values were higher. 

Finally, the presented methodology and the calibration results allowed us to find prices of European 
options currently listed on the EEX market. As the assumed 3-regime MRS model seems to be adequate to 
describe dynamics of electricity spot prices, the results of the paper can be used for reasonable pricing of 
electricity derivatives and, hence, yield an effective risk management tool. 
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Appendix 

Derivation of the option price formula. Using standard arguments the op tion price is the discounted ex - 
pected value of the payoff function under the pricing measure (see e.g. iMusiela and RutkowskiL Il997 ). 
Moreover, analogously to formula (TTUl) this expectation is equal to: 



C T {K) = e- rT E x [(B T - K) + \B ] = e - rT E x [(X T + g T - K) + \? ] = (32) 

= e- rT L^E x [(X T , b - K') + \F ] +pZ ] E [{Xy Tls K')+\F a ] +p$E [(X [TUd A") + |-F ] j, 

where K' = K — qt and (x) + = max(0, x). 

We start with pricing the base regime part Cx,b(K): 

/>oo 

Qr,b{K) = E x [(X Ttb - K')+\F ] = / (x - K')f XT bl x 0b (x)dx, (33) 

Jk> 

where fx Tb \x b (%) denotes the density of Xr.b conditional on Xqj,. From f|14|> and Ito's lemma we have 
that: 

X T .b = X^ b e-P T + % (1 - e^ T ) - f er^ T ^\{u)du + a b f e~^ T - u Uw x . (34) 



o 
Hence, Xt^ given Xq^ has a Gaussian distribution with mean 

and variance 



E(X T ^ b \To)=X ^ T + ^(l-e^ T )- I e-^ T -^\(u)du, (35) 



Var(X T , b \^) = ^(l-e-^ T ). (36) 



Denote the mean by m and the variance by s . We have: 

C T , b {K) = -i=- [°°(x - K') exp (- {X ~™ )2 ) 



2tts Jk> \ 2s : 

K! — m 



(K 1 - m)' , 
exp -^ — ' +{m-K') 



f2w V 2s 

14 



1-$ 



(37) 



Now, we turn to the pricing of the spike regime part. Observe that if A' < c s , then E[(X\ T \ S — 
if') + |J-o] = E(X\T\ a ) — A'. Assume that A' > c s . Denote the log-normal pdf by fLN(^ 3 ,a 2 ) an d the cdf 
by F LNiMs ^2y Then 

/>OC 

E[(X lTUs -K') + \T ] = / {x-K')f LN{ ^ al) {x-c s )dx = 

JK' 

"log(JJf # - c s ) - n„ - cr 21 



-r- y I <l-$ 



and we have that 



£[(X LrJ , s -A") + |.Fo] = V'>c}1 <*p hu 



(A'-c s )[l-F iJV(M8>CTf) (A'-c s )] 

log(A' - c„) - n a - a 2 s 



(38) 



1 -$ 



-(Jf / -c.)[l- J F iJv0l ., o2) (A"-c.)] 

^2 



■yl+C»- A'' 



+I{if'<c 3 } 

Similarly, we can price the drop regime part: 



E\(Xy Tld -K') + \Fo\ hK><c d }{ -exp(/i d + -^| 



(39) 



log(c d - A') - n d -aj 



(T,l 



+ (c d -K')F LN( ^ al) {c d -K') 



(40) 



Finally, letting C T ,i{K) = E[(Xy T \ A - K') + \jF ] for i £ {s,d} and combining formulas ((221), ([37]). (j3"9")l and 
(|4T)|) yields the result. 

Derivation of the forward price. First, note that 

E x {P T \F t ) = E\X T \F t ) +g T = E x (l {R[Tj=b} X T , b + I {R[T]=s} X lTl , s +I {R[Tj=d} X lTUd \T t ) + g T - (41) 

Since X\rj t i, for i = s, d, is independent of R[t] an d -T 7 *, we have: A A (I{ fl ^j-jXr^l-Ft) = P(R\t\ — 
i\J 7 t)E(X^ T j i ), i — s,d. Secondly, we have 

E x (l {RlTj=b} X Ttb \T t ) = E x (I {RlTj=b} \T t )E x (X T , b \T t ) = P(R m = b\F t )E x (X T , b \T t ) 
and from the base regime definition and Ito's lemma: 

E\X Tfi \T t ) = E\X t , h \Tt)e-K T -* + |(l - e~^ T ^) - £ 



-P(T-u) 



\{u)du. 



Moreover, from the definitions of the spike and drop regimes, see equations © and {7}, we get: 



E(X lTUs ) = e^ + ^s 



E(X lTUd )=c d -e^+^ 



and 



Finally, combining formula (gTj) with (|32|). (j43j) and (J44J) yields the result. 
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(42) 

(43) 
(44) 



Derivation of the price formula for an option written on a forward contract. We start the derivation of the 
option price formula with the following observation. 
If R[ t j = b, then the forward price is given by: 

X t . b I w{T 1 ,T 2 ,T)P{R VT] = b\R [ti = b)e-^ T - t) dT + 



f[Ti,T 2 ] 
Jt\{R w =b} 



■1\ 



f 2 w(T u T 2 ,T)P(R lTl =b\R\t\ =b) 



1 -e 



-P(T-t) 



-/3(T-«) 



X(u)du 



dT 



l-T-2 



+ ( e ^+K 2 +Cs ) f ' w (T u T 2 ,T)P(R lTi = s\R\t } = b)dT + (45) 

+(c-e w+ *°3) / 2 w(T u T 2 ,T)P(R [Tl = d\R [ti =b)dT+ f \(T x ,T 2 ,T)gT&r. 



Moreover, if k is such a number that Ru\ = i, -R|tJ-i 7^ ^> ■••) ^L*J- fc +! ^ b,R^ t ^_ k = b, for i £ {s,d} (i.e. 
the last base regime price before time t was observed fc-periods earlier), then the forward price is given by: 



/; 



[Ti,T 2 ] 
*|{-RL*J = ^- R L*J-i^'''--''- R LtJ-fc=''} 



X [t] _ k+hb / w(T 1 ,T 2 ,T)P(R lT] =b\R [t} = i)e-^ T -^ +k - 1 UT + 
* w(T u T 2 ,T)P(R m = b\R lti =i)x 



- (i - P -^ T -\*i+ k - i A - f 



s -/3(T-u) 



\{u)du 



dT 



T-2 



c s ) I w{T u T 2 ,T)P{R m =s\R\t\ =i)dT + 
+(c d -e^+^) f 2 w(T u T 2 ,T)P(R m = d\R [t] = i)dT + 

•M 



w(Ti,T 2 ,T)g T dT. 



(46) 



'i'i 



Formula ()45[) is a simple consequence of equation (|23fl and the fact that i?i t j =6 implies that E x {X t fi\Fi 



L*J 
X t> fc. In order to show (|4"6")l , observe that for .R^j = z, -Ritj— i] 7^ &; •••) -R[*J- fe = b, i = s,d, we have 

^(X t , 6 |J" t ) = ^[x LtJ _ fc+1 , 6 e-^-W+*- 1 ) + |(l-e-«*-W+*- 1 ))- 

e -«*- u >A(u)du + o* / e-"(*- u )dW„ A |X Lt j_ fc+1 , 6 

L*J-fe+i ^ L*J — fc+i 

= ^itJ-fe+i ^-^-L*^- 1 ) + £ (l - e -/3(*-L*J+fc-D) _ /' e -«*->A(tOdu. (47) 

Combining equations (|23[) and (|4"Tf yields the result. 

Now we can derive the option price formula. The option price is equal to the expected future payoff. 
Therefore, we have 



Cfl TuT2] (K) = e- rt E x \(fl Tl ' T2] -K) + \J'o 



(48) 



Observe that the forward price can be written as 



/, 



[Ti,T 2 ] 



L'J 



ATuT 2 ] 



T ati,t 2 ] sr^ st-t, f i -. 

"■{R M =b}J t \{ Rw= i,} I" 2^, Z^ !1 { R lt\=i,Rit\-i¥=b,---,Rit\-k=b}Jt\{R. vti =iM. vti _ 1 =tb,...,R lti _ k =b}i 
ie{s,d} fe=l 

16 



where /I^Lfti and f}J}£*} =iR: ^ ,^ h R _ , =hV , i € {.5,4, fc = 1,2,..., |_*J, are given by (g5]) and gSJ). 



Hence, 



/t|{R LtJ =6} 



t|{Rm=i,fluj_i#fe,...,-R| t j_A,=6} 



C/i Tl ' T2l (i^) 



-rtJ^A 



/ f [Ti,T 2 ] 
Wt|{R LtJ= fc} 



^) + |J- P(fl LtJ =b\R = b) + 



L'J 
I 

JS{s,d} fe=l 



Z Z {-^ [rfKi 1 



:-Rl*J-i^ 6 Ritj-k=b} 



K)+\T a 



:P(i? LtJ =i,i? LtJ _i ^6,...,iJ Lt j_ fc =6|i?o = &)}■ 



(49) 



Now, observe that 



E> 



if, 



[Ti,T 2 ] 
t\{R lti =b} 



K)+\JF 



(X tjb A (b) + B (b) - K) + | To 



A a {b)E } 



X, 



i.h 



K - B (b) 
Mb) 



T 



Mb)C t ,b 



K-B (b) 
Mb) 



(50) 



9t 



where Ct,b I A i° b \ + 9t ) is the 'base regime part' of the price of a European call option written on the 

electricity spot price with maturity t and strike — A ,° b \ + gt , see equation (|18p and A k , B k are defined in 
equations (|25 ]) -(|25 )) . Similarly, 



E* 



(ATi,T 2 ] _ v-\+\-p 



A k (i)E > 



X 



\t\-k+l,b 



K-B k ( i)\ + 
A k (i) 



l^o 



Mi) c Yt\-k+i,b [ — ■j-A K<?LtJ-fc+i ) 



for i G {s, d} Finally, combining formulas (|4"5j) . ([50]) and ([ST|) completes the proof. 
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